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Abstract 

We show that the point spectrum of the standard Coulomb-Dirac op- 
erator Ho is the limit of the point spectrum of the Dirac operator with 
anomalous magnetic moment H a as the anomaly parameter tends to 0. 
For negative angular momentum quantum number k, this holds for all 
Coulomb coupling constants c for which Ho has a distinguished self-adjoint 
realisation. For positive k, however, there are some exceptional values for 
c, and in general an index shift between the eigenvalues of Ho and the lim- 
its of eigenvalues of H a appears, accompanied with additional oscillations 
of the eigenfunctions of H a very close to the origin. 

1 Introduction 

In a 1930 letter to O. Klein ([Jjj letter 261) and in a survey to be given at the 
8th Solvay Congress, planned for October 1939, Pauli suggested to describe the 
motion of a particle with rest mass m > 0, charge e, spin h/2 and magnetic 
moment (l + a)/ig (hb = -fi- the Bohr magneton) in an electric field — V<& and 
a magnetic field B = curl A by means of the operator 

H = cat ■ Op— -A) + mc 2 [3 + e$ - a/i s (ia ■ V$ + a ■ B). 
c 

A revised version of this review appeared in |S] (the relevant equation is (91)); 
the original manuscript was not published until 1993 in ^Uj pp. 827-901. 

If $(x) = V(|x|) and A = 0, self-adjoint realisations of H in L 2 (R 3 ) 4 are 
unitarily equivalent to the orthogonal sum of self-adjoint realisations of 

mc / k a/jn ,, ,\ e 
H a = ^^ + {- + ^V'{r))a 1 + -V { r) 

in i 2 ((0, oo)) 2 , where k € Z \ {0} is the angular momentum quantum number. 

The first mathematical treatment of the operators H and H a is due to 
Behncke [2E1I3I- He showed that H a has a unique self-adjoint realisation if 
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a 7^ for a very large class of potentials V, including the Coulomb potential 
V(r) — for all values of the coupling constant Ze (for alternative proofs 
see jS] and PP). This is in marked contrast to the case a = where it is well 
known that H is essentially self-adjoint on its minimal domain if and only if 
(^r) 2 < k 2 — j- For larger values of Z the singular end-point is in the limit- 
circle case; but as long as (^r) 2 < k 2 one still has a distinguished self-adjoint 
realisation of H defined by the requirement that functions in the domain behave 
like the principal solution of the eigenvalue equation of H at 0. 

The location of the essential spectrum of H a is comparatively easy to de- 
termine and for a wide range of potentials, notably the Coulomb potential, 
coincides with that of H 0EH2. 

In the following, we normalise constants and write the Coulomb Hamiltonian 

as 

d (k a \ c 

H a = -«*2— + 03 + I - + -o ) O"! + ~, 

dr \r r z / r 

assuming a,c < 0. (The cases of positive c and/or a can be reduced to this 
situation by means of suitable unitary transformations.) The discrete spectrum 
of Ho accumulates at the right end-point of the gap (—1,1) in the essential 
spectrum. Since Cq°((0, oo)) 2 is a common core for Hq and H a if (and only if) 
c 2 < k 2 — j, H a converges to Hq in the strong resolvent sense as a — > ( jllj 
Thm VIII. 25(a)), and as a consequence, the spectrum of H a cannot expand in 
the limit a — > for this range of the parameters c, k. However, it is reasonable 
to expect (and has been used as a basis for a perturbative calculation of the 
eigenvalues of H a ) that the point spectrum is stable in the limit a — ► in the 
sense that the eigenvalues of H a converge to those of Ho, and each eigenvalue 
of Ho is the limit of exactly one eigenvalue branch of H a . Decoupling the 
eigenvalue equation of H a and using a comparison theorem for principal and non- 
principal solutions of second-order equations, Behncke 0] proved this stability 
for k 2 > c 2 + (|) 2 if K < 0, and k > c 2 + | if n > 0. He conjectured that 
k 2 > c 2 + | might be sufficient in the latter case. (Farther-reaching conjectures 
are to be found in P- 218 seq.) 

In the present paper, we study the convergence of the point spectrum of 
H a as a tends to 0, for the whole parameter range for which a distinguished 
realisation of Ho exists, i.e. for k 2 — c 2 > 0. We find a surprising qualitative 
difference in the limiting behaviour depending on the sign of k. Indeed, for 
negative k, the eigenvalues of Ho are exactly the limits of eigenvalues of H a for 
all values of c. For positive k, however, there are (finitely or infinitely many) 
exceptional values Co > c\ > . . . in (— k, 0); for c G (c m ,c m _i), the eigenvalues 
of Ho are still the limits of the eigenvalues of H a , but with a shift of size m in 
the eigenvalue numbers. This shift is reflected in the appearance of m additional 
oscillations of the corresponding eigenfunction of H a , compared to that of Ho, 
very close to the origin. It seems a delicate question to decide whether the 
number of exceptional values c m is finite or infinite; in any case it grows beyond 
all bounds with increasing k. 

More precisely, we have the following results. 
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Theorem 1.1 (Spectral convergence and stability for negative k) 

Let k < 0, c € (k, 0), and let Ao [not] be an eigenvalue of the Coulomb-Dirac 
Hamiltonian 

d k c 

Ho = —ia 2 — + cr 3 H fii + -. 

dr r r 

Let < e < dist(Ao, o-(H ) \ {Ao})/2. Then for a < with sufficiently small \a\ 
the Hamiltonian with anomalous magnetic moment 

d ( k a \ c 

H a = -«cr 2 — + 03 + - + — Ci + - 

dr \r r A / r 

has exactly one [no] eigenvalue X a in (Ao — e, Ao + e). 

Theorem 1.2 (Spectral convergence and stability for positive k) 

Let k > 0; then there are at least [(Klog4)/7r — 1] values > c > c\ > . . . > — k, 
which can only accumulate at —k, such that the following holds. 

Let c S (— n, 0) \ {cq, ci, . . .}, and let Xq [not] be an eigenvalue of H . Let 
< e < dist(Ao, cr(i?o) \ {Ao})/2. Then for a < wit/i sufficiently small \a\, H a 
has exactly one [no] eigenvalue X a in (Xq — e, Ao + e). 

Here [x] := sup{m £ Z | to < a;} denotes the Gaufi bracket. 

The proof of these theorems is based on oscillation theory, in particular on 
an asymptotic study of the behaviour of the Priifer angle of solutions of the 
eigenvalue equation for H a (i.e. the solutions of equation J2J below) as a tends 
to 0. After rescaling g = r/\a\, the mass term and spectral parameter are lower 
order terms in the limit and can be omitted in order to obtain an overview of 
solutions in the asymptotic regime near the origin. 

The direction field of the resulting simplified equation Q (where k = 
and a — — sgnK), and hence the qualitative behaviour of its solutions, shows a 
fundamental difference depending on the sign of a. 

For negative k (a = 1), the (g, $) plane is divided into essentially horizontal 
strips in which the right-hand side of Q is alternatingly positive and negative 
(cf. Fig. as a consequence, the distinguished angle i9o which corresponds to 
an £ 2 (0, •) solution of the eigenvalue equation for H a , cannot change by more 
than 7r, and eventually tends towards an asymptotic value "$+(c) (Prop. 2.3 a) 
which turns out to be the limiting angle at of the principal solution of the 
eigenvalue equation for Ho as well (Lemma 3.1 b). On the original r scale, the 
convergence becomes faster as a — > 0. 

A stability argument (Lemma 3.2, Prop. 3.3) then shows that for a certain 
point R > the influence of the previously neglected mass and spectral param- 
eter terms can be controlled on (0, R), and that the solution of the full Priifer 
equation J2J) converges at that point to the Priifer angle of the principal solution 
of the equation for a = 0. 

Theorem 1.1 follows in view of the uniformity of this convergence with re- 
spect to the spectral parameter, and the fact that the presence of the anomalous 
magnetic moment does not significantly affect the behaviour of the solutions at 
00. 
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Figure 1: Zones in the directional field of the simplified equation for a — 1. The 
arrows indicate the sign of the right-hand side of QJ, the solid lines represent 
its zeros. 



A curious phenomenon occurs, however, in the case of positive k (a = — 1). 
For q close to and for large g, one again has ^-regions of opposite sign of the 
right-hand side of JJJ , and hence of essential confinement of the solutions, and 
for large g, the distinguished angle i?o generically tends to the limit $+(c) mod tt. 
In contrast to the previous situation, there is now a p-interval (g_(c), g+{c)) on 
which the right-hand side of (JTJ has no zeros and is strictly negative (cf. Fig. 
0). 

Depending on c, the size of this gap increases, vanishing as c — ► and 
becoming infinite as c — > —k. Moreover, one can show that the angle $o changes 
by several multiples of tt in this interval if — c is large enough (Prop. 2.3 b), 
and hence will eventually converge to $ + (c) — rmr for some m € No- Thus the 
corresponding £ 2 (0, •) solution of the eigenvalue equation for H a will perform m 
oscillations which are absent in the principal solution of the eigenvalue equation 
for Hq, to which it however converges in phase at the point R > 0. 

Theorem 1.2 then follows by way of the same stability argument as before 
(Lemma 3.2, Prop. 3.3). 

As c crosses the exceptional value, a transition of the asymptotic limit from 
j?_l_(c)— rn7r to (c) — (m+l)7r takes place. At c = c m , the distinguished solution 
i?o of the simplified equation Q approaches an unstable limit as g — > oo. The 
asymptotics of the solutions of the full Prufer equation appear to be a rather 
delicate matter in these cases, and the limiting behaviour of H a with c = c m 
remains an interesting open question. 

The basic analytical tool of the present paper is the study of the direction 
field of Prufer and Riccati type ordinary differential equations near a singularity 
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Figure 2: Zones in the directional field of the simplified equation for a = — 1. 
The arrows indicate the sign of the right-hand side of (JTJ , the solid lines represent 
its zeros. 



(for related earlier but simpler results see [Jj and the references therein). The 
underlying comparison techniques are outlined in the Appendix. 



2 The simplified equation 

In this section we consider the scaled and simplified equation 

Q-d' = c+ (k + a/g) sin 2^ (1) 

(with k > 0, c E (—k, 0) and a E {—1, 1}) which arises from the Priifer equation 
(J5J equivalent to the eigenvalue equation for H a by omitting the 0(1) (g — > oo) 
terms after rescaling r = \a\g, which eliminates \a\. 

The asymptotic zeros r d± (c) of the right-hand side of equation satisfy 

< i?_(c) < tt/4 < <& + (c) < it/2, 



sin2^±(c) = -c/k, tani9±(c) = MVMES. 

Moreover, in case a = — 1 there are the exceptional points g±(c) = l/(k ± c), 
between which the right-hand side of QJ is strictly negative; for later convenience 
we define g+(c) = if a = 1. 

We first show that has a stable and an unstable asymptotic critical point 
at both singular end-points and oo (Lemma 2.1, 2.2). The distinguished 
(unstable) solution at 0, i?q , converges to the stable limit at oo for all c if a = 1; 
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if a = — 1, one has the convergence with a shift by an integer multiple of ir 
unless c is one of a sequence of exceptional values (Prop. 2.3). 

Lemma 2.1 (The existence of a distinguished solution at oo) 

For each c € (— k, 0), 0) has a unique solution i?oo('jc) such that 

lim i9oo(e, c) = i?-(c). 

£ — >oo 

oi/ier solutions $(-,c) satisfy either i)(g,c) = #oo(f?, c) + m7r (g > 0) 
/or some to G Z, or eZse lim^oo $(g, c) = $+(c)mod7r. Furthermore, for each 
g > 0, #oo(f?> ') w continuous and strictly decreasing. 

The proof of this lemma, based on an asymptotic study of the direction field of 
l|T[l. can be found in the Appendix. Similarly, one can prove 

Lemma 2.2 (The existence of a distinguished solution at 0) 

For each c £ (— k, 0) there is a unique solution i9o (■, c) such that 

i- n / \ Ft if oi = — 1, 
limtfoG?,c) = j/a = L 

j4/Z other solutions are either shifts of *&q by an integer multiple of n, or have 

,. a/ \ |7r/2mod7r ifa — —l, 
e ^o v ' [0mod7r ifa = l. 

For fixed g > 0, "&o{g, •) *s continuous non- decreasing. 
Proposition 2.3 

a) If a — 1, we have lim $o(g> c) = $+(c) (c € (— k, 0)). 

b) If a = —1, there are at least [(fclog4)/7r— 1] values > cq > c\ > C2 > . . . > 
—k (accumulating at —k if infinitely many) such that 

lim t?o(g) c) = i?+(c) - m7r (c € C m ), 
and lim r do(g,c m ) = i?_(c m ) — mix. 

Q — >00 

Here C m := (c m , c m _i) (where c_i := 0J /or aZZ to S No /or which c m exists; 
if there is a minimal c mmax > —As, define Cm max +i := (— k, c rn . msa ). 

Remark. For later convenience, we also define Co := (—As, 0) in the case a = 1. 

Proof, a) The interval [7r/4, 7r] is stable for by Lemma 4.2 (cf. the Appendix), 
so i9o (■, c) (which by Lemma 2.2 is close to ir/2 for small g) cannot tend to 
i9_(c)mod7r. The assertion follows by Lemma 2.1. 
b) In the limit c — > 0, we have, using fl'(g) > (l/g)(c— \k — l/f?|) 5 

i?o(fl+(c), c) - o (fi- (c), c) > § (c - i + A;) dg + J$ c) \{c+\-k)dg 

= k ~^Tc)-( c + k ) lo §(^- (c)) - + fc + (c - fc) log(fco+(c)) 
= ( c -fc)log^-(c + A ; )log^^0. 
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As a consequence, lim i?o(f?> c) = #+(c) for c < sufficiently close to 0. 

On the other hand, noting that 'd'(g) < (f /g)(c+\k— l/g\) and $ OCl (g + (c), c) e 
[0, 7r/4], we have 

0o(l/fc, c) - t?oo(l/A!, c) < i?o(ff-(c), c) + J e V ( c) i (c + i - k) dg 
-$ 00 ( e+ ( c ),c)+ j$ c) \{c-\ + k)dg 
^ ^ - * + T^J + (* ~ c ) l0 S( fc ^- ( c )) + 7^j- k +( k + c ) ^g(kg+(c)) 

7T — fclog4 (c — > — fc). 

The assertion follows, as ^(l/k, •) and ^(l/fc, •) are continuous and monotone 
non-decreasing and decreasing, resp. □ 

Remark. The first part of the proof of Proposition 2.3 b) yields a quantitative 
estimate for the first exceptional value cq. Indeed, using Lemma 4.1 one can see 
that i?o(£-(co),co) > x and i? (£+(co), c ) < f , so 

77 < (c - fc) log fc + c ° - (c + fc) log C ° . 
2 k k 

Setting x := —co/k e (0, 1) and observing that 

xiogi+f -bga-x 2 ) =Y ,«L lX v + 

< t-2 fo, x 2 7\ _ x 2 (18~llx 2 ) 
— X \° ~r T^x 1 6 J ~ 6(1-k 2 ) ' 

we find llx 4 - (18 + 3ir/k)x 2 + 3ir/k < 0, and thus 

Co > (18fc 2 + 37rfc - v/(I8fc 2 + 37rfc) 2 - 1327rfc 3 ). 

For fc = 1, this gives the bound c < —0.64157. 

Hence, one has convergence and stability of the eigenvalues by Theorem 1.2 
for all k € Z \ {0} at least for nuclear charge number Z < 87. 

3 The convergence of the original equation for 

a — > 0, c 6 C m 

Throughout this section, we fix a,fi e { — 1,1} and fc > 0, and assume that 
ceC m for some admissible m e Nq. 
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Consider the Priifer equation 



6' = - + (- + ^)sin2e + /icos2e-A (2) 
and the corresponding equation for a = 0, 

X' = - + - sin 2X + /i cos 2X - A. (3) 

r r 

We shall study the following distinguished solutions, whose existence and 
properties can be proved in analogy to Lemma 2.1. 

Lemma 3.1 

a) There is exactly one solution 0o of 0) with 

i- ^ i \ \ f n "if a = —1 

hme (r,c, a ,A) = | 7r/2 . fa=1> 

all other solutions either being shifts ofOo by an integer multiple of it, or having 

_j7r/2mod7r if a = — 1 



lim 9(r, c, a, A) 



mod 7r if a = 1. 



For fixed r, c and a, &o(r, c, a, •) is continuous decreasing, 
b) There is exactly one solution Xq of |3J) with 

limX (r, c, A) = #+(c), 

all other solutions either being shifts of Xq by an integer multiple of it, or having 

lim X(r, c, A) = (c) mod n. 

For fixed r and c, X (r, c, •) is continuous decreasing. 

Now we show that the solutions X and 0o become asymptotically close mod tt 
at some point for small a; this is a consquence of the convergence of the solution 
$o (which is close to Oo) of the simplified equation (Prop. 2.3). 

Lemma 3.2 For each e > there are ro(e) € (0,e] and ao > such that for all 
A G [—1, 1] and — ao < a < 

|6o(r (e), c, a, A) + run - tf + (c)| < f , 

|JC (n,(e),c,A)-i? + (c)| <§, 

and consequently |0o(ro(e), c, a, A) + m7r — Xo(ro(e), c, A)| < e. 

Proof. Let 7 > be so small that [0 — 7,0 + 7] C C m . Then there is S £ (0, 7) 
such that \d + (c±5) — #+(c)| < e/3. By Proposition 2.3 there is go > such that 
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|#o(f?> c± 5) — fl + (c± 5) + nm\ < e/3 (g > go), so the function r h- >• ■do(r/\a\,c) 
(which is a solution of the simplified equation 

A// n c ,k a\a\ , . 
9'(r) = - + (- + -^j 1 ) sin 26) 

satisfies |i? (r/|o|,c± 5)- t? + (c) + m7r| < 2e/3 (|o| < > 0). Now let 
r (e) := min{(5/2, e}. Estimating in (j3Jl, @ 

|/ucos2X - A|, |^ cos 29 - A| < 2 < - (r < r (e)), 

r 

Lemma 4.1 yields the bounds 

tf (r/|a|, c - S) < 6 (r, c, a) < tf (r/|a|, c + 5) (0 < r < r (e)) 

and 

i?+(c-<S) <X (r,c) <tf + ( c + £) (0<r < r (e)), 

and the assertion follows. □ 

The preceding lemma shows that the solutions converge to each other; however, 
the point of comparison r (e) depends on e and tends to rather rapidly. We 
now show that the convergence remains stable, and hence also holds at a certain 
fixed point R. 

Proposition 3.3 There is R > such that lim Qo(R, c, a, A) = Xq(R, c, A) 

a— vO 

uniformly w.r.t. A £ [—1,1]. 

Proof. Consider the Riccati equations for z(r) := tanATo(r, c) and y(r,a) 
6 (r, c, a), 

r / = (c - (A + l)r)z 2 + 2kz + c - (A - l)r (4) 

and 

ry' = (c- (A + l)r)z 2 + 2{k + a\a\/r)z + c - (A - l)r. (5) 



Let j/+(c) := tantf+(c) = . Choose < d < y/k 2 - c 2 /(3\c\) and 

— c 

< R < d so small that 

3 ||A + l|(y+(c) + d) 2 + |A- 1|| i2,6|A + l|(j/+(c) + d)R < Vk 2 - c 2 
for all A S [— 1, 1]. For any f S (0, i?), set 



/ », Vk 2 — c 2 A 

ai(r) := p-r r. 

6(2^1 + 1) 

Now let e > 0, e < R < d. By Lemma 3.2, there is ro(e) < e and ao > 
such that 

|«(r (e)) - y+(c)| < e/2, |j/(r (e), a) - y+(c)| < e/2 (0 < |o| < oo). 
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On [ro(e),R], the interval [y+(c) — d,y+(c) + d] is stable for Q and |J5J by 
Lemma 4.2 if |a| < ai(ro(e)), since the right-hand side for y+{c) ± d takes the 
value 

(c — (A + l)r)(2/+(c) ± d) 2 + 2(fc + a\a\/r)(y+(c) ±d) +c — (X — l)r 
= cd 2 ± 2(k + y+(c)c)d + 2^(y+(c) ± d) - (A + l)r(y+(c) ± d) 2 - (A - l)r 
Td \2VW^ ±cd± 2 2M(i±M ± i) T ((A + l)(y + (c) ± d) 2 - A - , 



and the factor in square brackets is not less than \Jk 2 — c 2 . 

Hence |y(r, a) - 2/+(c)|, |z(r) - y+{c)\ < d for all r <E [r (e),i?] if < |a| < 
min{a ,ai(r (£))}. 

In the differential equation for the difference x(r, a) := y(r, a) — z(r), 

rx' = (c- (A + l)r)(y + z) + 2fc) .x - (6) 

r 

the factor of x on the right-hand side can be estimated 
c(y{r,a) + z{r)) +2k - (A + l)r(y + z) 



< -2(Vk 2 -c 2 + cd) + 2R\X + l|(y+(c) + d) 



< -Vk 2 - c 2 . 



Hence, if a additionally satisfies 



Vfc 2 - c 2 r (e) 
1 1 2(y + (c)+d) 

then the interval [—e, e] is stable for © on [ro(e), R] by Lemma 4.2, and because 
of \x(ro(e), a)\ < e we find |y(i?, a) — -z(i?)| < e. □ 

For the proof of the main theorems, we need the following distinguished solution 
of 10 at infinity, whose existence and properties can be proven along the lines 
of Lemma 2.1. Let X±(\) £ [0,tt] such that cos 2X± (A) = fiX and sin 2X± (A) = 
±nVT^~X 2 (A £ [-1, 1]). 

Lemma 3.4 For X £ [—1,1], i/iere is exactly one solution X^ of J3J) with 

lim X 00 (r,c,A) = X_(A), 

r — >oo 

all other solutions either being shifts of Xq by an integer multiple of it, or having 
lim X(r, c, A) ~ X + (X) modn. 

r — >oo 

For fixed r and c, X^ (r, c, •) is continuous increasing. 
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Now we are in a position to prove the main results. 

Proof of Theorem 1.1. 

Let a = 1, fi = —1, k = \n\ = —k, and R > as in Proposition 3.3. It is 
sufficient to prove the assertion for the auxiliary Hamiltonian 



d (n a \ c 

H a = + 0" 3 + I V -jX(0,R)(r) ) Ol + - 

ar \ r r / r 



instead of H a , since the eigenvalues of H a are within |a|/i? of those of H a . 
Introducing the Priifer transformation 

, , / sin 

U = \U\ Q 

1 1 \_ — COS V 

in the eigenvalue equation (H a — A) u = 0, we find the Priifer equation for $ 

i9' = - cos 2tf - (~ + ^X(o,h) (r)) sin 2d + °- - A, 

which in view of the above choices for a, /i and k coincides with on (0, R) 
and with J2J on (_R, oo). 

For the Priifer radius \u\ we have 

(log|u|)'(r) = - sin 2d+ + ^X(o,Ji)(r)) cos 2,9. 

Hence, if § = Oo on (0, R), where Oo is the distinguished solution from Lemma 
3.1 a) (with a = 1), we find (log|ii|)'(r) ~ |a|/?" 2 (r — > 0), and hence 



\u\ 2 (r)dr ~ const / e~ 2|a|/r dr< 



oo. 



Thus Go is the Priifer angle of an £2(0, •) solution of the eigenvalue equation 
for H a . 

Similarly, for a = we find for a solution u with Priifer angle d — Xq, where 
Xq is the distinguished solution from Lemma 3.1 b), 



(logM)'M 



Vfc 2 - c 2 



so \u\(r) ~ const r^ - ^ (r — > 0), whereas for all other solutions v 



(logM)'(r) ~ 



and hence \v\(r) ~ const r~^ k2 ~ c2 (r — > 0). Thus X is the Priifer angle of the 
principal solution of the eigenvalue equation for Hq. 

Analogously, if (for either a < or a = 0) — Xoo on (R, 00), where Xao is 
the distinguished solution from Lemma we have 

k 

(logld)'(r) = -31112X00 cos2Xoo VI - A 2 (r 00), 

r 
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so 



/ \u\ 2 (r)dr ~ const / e^ 2y/TI ^ r dr < oo. 

Thus Xoo is the Priifer angle of an L 2 (-, oo) solution of the eigenvalue equations 
for H and H a . 

As a consequence, the eigenvalues of H a are the (isolated) values of A at 
which the monotone decreasing continuous function Q$(R, a, c, •) and the mono- 
tone increasing continuous function X oc (i?, c, •) take the same value mod 7r. 
Similarly, the eigenvalues of Hq are the intersection points mod7r of Xq(R, c, •) 
and Xoo(i?, c, •). Hence the assertion follows in view of the uniform convergence 
of Oo(i?, a, c, •) to X (R, c, •) as a — ► (Prop. 3.3). □ 

Proof of Theorem 1.2. Let a = —1, /i = 1, k = \k\ = k, and R > as in 
Proposition 3.3. As in the preceding proof, it is sufficient to show the assertion 
for the auxiliary Hamiltonian H a . We now use the Priifer transformation 

i | ( costf 
u = kt . 

which leads to the Priifer equation 

i?' = cos 2d + (~ + ^X (0 ,fl) (r)) sin 2d + ~ - A, 

which in view of the above choices for a, [i and k coincides with |[5J on (0, R) and 
with J3J on (i?, oo). By studying the asymptotics of |u| as above, we again find 
that Oo from Lemma 3.1 a) corresponds to an L 2 (0, •) solution of the eigenvalue 
equation for H a , Xq from Lemma 3.1 b) corresponds to the principal solution 
of the eigenvalue equation for H$, and Xoo corresponds to an L 2 {-, oo) solution 
of cither eigenvalue equation. 

Hence the assertion follows as in the preceding proof. □ 



4 Appendix 

In the proofs of this paper we frequently use the following fundamental obser- 
vations about first-order ordinary differential equations. The first lemma (cf. 
[H] p. 27) is sometimes called Caplygin's inequality, but actually goes back to 
Peano (0 p. 44). 

Lemma 4.1 Let I C K be an interval, x$ £ I and fj : I x K — > K locally 
integrable in the first, and locally Lipschitz continuous in the second argument, 
j € {1,2}, with f\(x,y) < f2{x,y) {x G I,y S R). Furthermore, let y® < y®, 
and yj be the solution of the initial value problem 

y'(x) = f j (x,y), y(x ) = y° (je{l,2}). 

Then yi(x) < y2{%) {x £ I,x > xq). 

An immediate consequence is the following stability criterion. 
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Lemma 4.2 Let I C ^ be an interval, f : I X R — ► K locally integrable in 
the first, and locally Lipschitz continuous in the second argument. The interval 
[yi > 2/2] is stable for the differential equation 

y'(x) = f(x,y) (7) 

ok I iff(x,yx) > 0, f(x,y 2 ) < (x G I). 

Here an interval J is called stable on J for J7J) if y(xo) G J y(x) e J 
(a; G /, x > xo) for all xq £ /. 

We now use these observations to prove Lemma 2.1. 

Proof of Lemma 2.1. 

a) Let cG (-fc,0). Define J± := {§ G M | ±(c + fcsin2i?) > 0}, and := {d G 
J±\W- #±(c)+jTr\ > e(j G Z)} (e > 0). 

Then for each e > there are 7 > 0, Po > g+(c) such that for g> Po 

±(c+ (k + a/g)wx'&) > 7 (i9 G J±^)- 

Now consider a solution #(•, c) of QJ. If for all e > there is P > and 
j'eZ such that 

|0( ft c)-0_(c)|<e (f?>P), 
this means that lim $(g, c) = Q-{c) — jn for some j G Z. 

£ — >oo 

Otherwise, let e be already so small that this is not true, and 7, Po > as 
above. Then there is go > Pq such that {^(go, c) — J?-(c) +J7r| > e (j G Z), and 
hence 

#(f?o, c) G j| e) U Ji £) U (J [0+(c) - jtt - e, i?+(c) - jtt + e]. 
In we have ±"&'(g, c) > 7/0, so if "&(go, c) G J± , then 

±i?(e,c) > ±i9(g ,c) - 7 log £>/ go 

as long as c) remains in J± ; consequently there is g\ and j G Z such that 
i?(£>i, c) G [i?+(c) - J7T - e, i? + (e) - + e]. 

The latter interval is stable by Lemma 4.2. As e > was arbitrary, it follows 
that lim i)(g, c) = i?+(c) mod7r. 

b) Let r be an open interval with T C (— fc.O), and g > g+{c) (c G T). For 
1) £ I and c G T, denote by -d(g, c, %9) the solution of (JTJ with initial value 

c, #) = 7?, and consider the sets 

Si := {{cj) G (-fc,0) x R I lim <&(g,c,3) = tf+(c) - jV}, 

^ — >oo 

j G Z. is open. Indeed, let (co,i?o) G Sj, S := ($+(c - 7r/4)/2. Since i?+(c) 
depends continuously on c, there is e > such that 

[0+(c) - <5, i?+(e) + «J] c (tt/4, tt) (c g [c - e, c + e]). 
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Now let go > Q be so large that | -i? ( g> , c o, $o) — $+(co) < 5/2. As d{g ,c, i>) 
depends continuously on (c, there is positive i < e with 

|i?(£o,c,i?)-i? + (c )+j7ii <5 (ce [c -£,c + £],i?e [i9 - £, + <?]), 

and therefore 1}(qo, c, 3) + jir lies in the stable interval (n/A, n). 
By a) lim d(g, c, i?) = i?+(c) — j7r, i.e. 

(c, i?) e Sj (c e [c - i, c + e], i? e [i? - e, + £])■ 

c) For each c e (—A;, 0), there is exactly one solution i?oo( - , c) such that lim i?oo(£i c) 

#-(c). Indeed, there is at least one, as c is in some suitable set T and the cor- 
responding sets 

s j (c) = {tfe(0,7r/4)| M)e^-}, 

J G {0, 1}, arc nonempty by a) and open by b). 

Assume there are two solutions #i < i? 2 with lim #j(f>) = $-(c), i G {1, 2}. 

Then </? := 2(tf 2 - > satisfies 

Qip'(g) =2(/c + a/ £ »)(sin2^ 2 -sin2?9i) 

= 2(k + a/g)(sm 2^ £2£ipl + cos 2i?j p 



- 2V/c 2 - c 2 v? > (g — ► oo), 
contradicting 93(g) — > 0. 

d) Finally, for fixed g > 0, $oo(£?>c) is strictly monotone decreasing: if c\ < c 2 
and i?oo(^, Ci) < floo{g, c 2 ), then by Lemma 4.1 

i?_(ci) = lim i?oo(eiCi) < lim ^(g, c 2 ) = i?-(c 2 ), 
which is not true. 

Also, $oo(g,') is continuous. Indeed, if c n /* c (n —* 00), then $oo(f?, c n ) > 
$oo(g,c) because of the monotonicity, so lim $00(0, c n ) > i9oo(g,c). However, 

n^oo 

'>' would imply (c, lim i?oo(g, c„)) € So in contradiction to the facts that So 

is open and (c„, ^(g, c„)) ^ S (n e N). The right continuity follows in the 
same way. □ 

Acknowledgement. H.K. wishes to thank E.B. Davies for his hospitality at King's 
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